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The linear response to a space- and time-dependent external disturbance of a system of dilute 
condensed composite bosons at zero temperature, as obtained from the linearized version of the time- 
dependent Gross-Pitaevskii equation, is shown to result also from the strong-coupling limit of the 
time-dependent BCS (or broken-symmetry RPA) approximation for the constituent fermions subject 
to the same external disturbance. In this way, it is possible to connect excited-state properties of 
the bosonic and fermionic systems by placing the Gross-Pitaevskii equation in perspective with the 
corresponding fermionic approximations. 
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Ultracold atomic Fermi gases are currently being inten- 
sively studied experimentally, as their low-temperature 
properties can shed light on fundamental questions re- 
garding degenerate Fermi systems. Froducing a paired 
superfluid is a particular challenge, and is being pursued 
by several groups^"'^. Since the attractive interaction re- 
sponsible for superfluidity in a Fermi gas can be tuned to 
reach the strong-coupling condition, it also appears possi- 
ble to study experimentally the evolution from superfluid 
fermions to condensed composite bosons. (These form as 
bound-fermion pairs in the strong-coupling limit of the 
fermionic attraction.) This would permit for the first 
time a controlled crossover between two basic quantum 
systems sharing the same spontaneous broken-symmetry 
(superfluid) behavior. 

Dilute condensed bosons have already been studied ex- 
tensively both experimentally^ and theoretically^. The 
diluteness condition is matched in almost all current ex- 
periments. Correspondingly, the Gross-Fitaevskii (GF) 
equations^''' have proven sufficient to describe the ob- 
served phenomena. This is true both for ground-state 
properties (via the time-independent GF equation) and 
for the dynamical response of the condensate to an ex- 
ternal disturbance (via the time-dependent GF equation 
in its linearized form®). 

With the aim of connecting properties of superfluid 
fermions and condensed composite bosons, the time- 
independent GF equation for the condensate wave func- 
tion was obtained in a previous paper^ as the strong- 
coupling limit of the Bogoliubov-de Gennes equations^° 
for superfluid fermions at zero temperature. 

In this paper, a corresponding connection is estab- 
lished between the linearized form of the time-dependent 
GF equation (from which the dynamics of the bosonic 
condensate can be obtained) and the strong-coupling 
limit of the linear response for superfluid fermions with 
a BCS ground state, subject to the same external dis- 
turbance. We will explicitly show that the fermionic 



density and current correlation functions map, in the 
strong-coupling limit of the fermionic attraction, onto the 
linear-response results for condensed composite bosons. 
The fermionic correlation functions are obtained from 
the time-dependent BCS (or broken-symmetry RFA) ap- 
proximation, while the linear-response results for the 
composite bosons are obtained from the time-dependent 
GF equation. Thus a connection between ground- 
and excited-state properties of the condensate both 
in the weak- and strong-coupling limits is established. 
On physical grounds the connection between the time- 
independent and time-dependent versions of the GF 
equation for condensed composite bosons and the BCS 
and time-dependent BCS approximations for superfluid 
fermions, rests on the zero-point motion in a quantum 
theory implying an oscillatory spectrum. The imprint 
of the quasi-particle spectrum should be found in the 
ground-state wave function^. 

The response functions of a system of condensed 
bosons described by the wave function <I>(r, t) can be ob- 
tained from the time-dependent GF equation^''': 
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in the presence of space- and time-dependent vector Acxt 
and scalar Vcxt external potentials. Here, c is the ve- 
locity of light, rriB and qs are the mass and charge of 
the bosons, and C/o = '^1^0,3 /ruB is the short-range re- 
pulsive potential expressed in terms of the bosonic scat- 
tering length aB {fi ~ 1 throughout). Note that we have 
generalized the time-dependent GF equation to include a 
vector potential, as a mere formal device to obtain, in ad- 
dition to the density-density correlation function already 
considered for neutral bosons,® the current-density and 
current-current correlation functions. In this way we can 
compare these functions with the corresponding quan- 
tities for fermions obtained within the time-dependent 
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BCS approximation in the strong-coupling limit. Accord- 
ingly, no Coulomb repulsion between bosons is included. 

The linear-response solution of Eq. (1) for Aext = 
is obtained^ by taking V;xt(r,t) = VqeM'"-'^*'' + 
V* e-'(qT-i^«) wave vector q and frequency uj, and 
by setting to linear order 



$(r, t) = [1 + u^{r)e-'^* - w„(r)*e 



(2) 



riB is the ground-state bosonic density, hb = UqUb the 
bosonic chemical potential, and the "small" components 

Ucj{r) and Vu;{r) are correspondingly taken of the form 



Wu(r) 



e*'! ''. The time- 



dependent GP equation (1) then reduces to a linear sys- 
tem for the unknowns Uq^,^ and Vq^,^. The solution is: 



*q,aj 



2 ^) ''^q,w — 



oj^-EBidy 



Vq. (3) 



q^/(2mB) is the free-boson energy and 



EB{q) = ^yE%{q)^ + 2^iBE%{q) the Bogoliubov disper- 
sion relation. The induced particle density and current 
are then obtained by expanding the expressions n(r, t) = 
|$(r,f)|2 andj(r,f) = Im[$(r, t)*V$(r, i))/mB] to linear 
order in Vext- One obtains: 



Sn{r,t) = 
j(r,i) = 



— 2 F / ^2 ^ext(r,^) (4) 

TUB - EBiqY 



riB 



.2 -ETJ-^VeM^^t) 



TUB - Esi^y 



(5) 



These satisfy the continuity equation dSn{r,t)/dt + V • 
j(r, t) = 0. By linear-response theory, the density-density 
and current- density correlation functions can be identi- 
fied: 



Xnn(q,w) 



riB 



Xj„(q,w) = 



ms ^2 - EBiq)"^ 
qw 



(6) 
(7) 



with the typical Bogoliubov form. 

Similarly, the linear-response solution of Eq. (1) 
when Vcxt = is obtained by taking Aext(r,i) = 
Aq e'(i '"-'^*) + A* e-'fi "-"*) and using the same form 
(2) for the wave function. Working in a gauge where 
V • Aext(r,t) is nonzero (such that q • Aq ^ 0) and re- 
tracing the above arguments, one obtains 



E%{ci) + 2nB+i^ Qb 



E%{q)+2nB-io QB 



Eb{ci) 



2mBC 



q-^q 



(8) 
(9) 



in the place of (3). The ensuing induced density and 
current now take the form: 



5n{r, t) = 



riB 



Loq 



ruB - EB{qy 
2nB E%{q)+2fiB 



(-^Aext(r,i)) 



(2ms)2 



^[q)+2HB (_qB_ 
-Esiqr^'y c 



Aext(r, t) 



(10) 
• (11) 



These replace (4) and (5), with the continuity equa- 
tion being satisfied once the diamagnetic contribution 
— ('J'B'^B)/(™Bc)Aext(r, is added to the induced cur- 
rent (11). The density-current and current-current corre- 
lation functions also acquire the typical Bogoliubov form: 



Xnj{q,oj) 

Xjj(q,w) 



riB 



Loq 



niB ^ 



EB{qf 



UB HB + E%{q)/2 



q ■ 



(12) 
(13) 



The correlation functions (6)-(7) and (12)-(13) describe 
the linear response of the bosonic condensate to an ex- 
ternal disturbance in terms of Bogoliubov quasi-particle 
excitations. 

We now show that the same form of the correlation 
functions is obtained from the fermionic time-dependent 
BCS (or broken-symmetry RPA) approximation in the 
strong-coupling limit of the fermionic attraction, when 
composite bosons with ijib = 2m and qB = 2q form. 
m and q are the fermionic mass and charge. We first 
express the fermionic correlation functions in the broken- 
symmetry phase in terms of the two-particle correlation 
function 

L{1, 2; 1', 2') = (T, [V'(l)V'(2)^t(2')^t(i')]) 

-g{l,l')g{2,2') . (14) 

Here, (• • •) is a thermal average, Tt the time-order op- 
erator for imaginary time t, 1 = (ri,ri,£i) with the in- 
dex £ = (J, II) identifying the components of the field 
operator in the Nambu representation (V'/(r) = V't('^) 
and Vi7(r) = V|(r)), and g{l,2) = - (T, [^(1)^^(2)]) 
is the single-particle Green's function. The two-particle 
correlation function (14) can be expressed in terms of the 
many-particle T-matrix^^: 

- L{1, 2; 1', 2') = g{l, 2') gi2, 1') + J d3456 

xe(l,3)e(6,l')T(3,5;6,4)e(4,2')e(2,5) . (15) 

An explicit form for T results from the Bcthe-Salpeter 
equation once a choice for its kernel is made. 

Within the (off-diagonal) BCS approximation we are 
considering, the Bctho-Salpotcr equation for T can be ex- 
plicitly solved in the frequency and wave- vector represen- 
tation. For fermions interacting via an attractive short- 
range potential regularized in terms of the two-body scat- 
tering length ap (with as = 20^)^^, the only nonva- 
nishing elements of the many-particle T-matrix corre- 
spond to the Nambu indices reported in Fig. 1(a). In 
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the strong-coupling limit (/3/i ^ — oo with inverse tem- 
perature (3 and fermionic chemical potential /i) of the at- 
tractive fermionic interaction, one obtains the following 
expressions for the nonvanishing elements of the many- 
particle T-matrix: 



rii(q,fi.) =T22(q,-f7,.) 

^ Stt iVL^ + / {2mB) jJ-B 

Ti2(q,0.) =T2i(q,f2,) 

^ Stt 



(16) 
(17) 



ill, = 2i/7r//3 {i^ integer) is a bosonic Matsubara fre- 
quency, fj,B = 2^ -\- is the bosonic chemical potential 
(with Eq = {ma?p)~^ the two-body binding energy), and 
Eb(s\) coincides formally with the Bogoliubov dispersion 
of Eq. (3). The fermionic single-particle Green's func- 
tions are correspondingly given by the BCS expressions 



C/ii(k,w„) = - ^?22(k, -w„) 
^2i(k,a;„) = ^i2(k,w„) = 



4(k) + ii^n 

A 



Ef^(k)2+w2 



(18) 
(19) 



ujn — (2n -I- {n integer) is a fermionic Mat- 

subara frequency, ^(k) = k2/(2m) — fi, and -E'(k) = 
-\/^(k)2 + for an (s-wave) isotropic gap function A. 

I --—I ' n 11^ 
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FIG. 1. (a) Nambu indices associated with tlie four nonva- 
nishing elements of the many-particle T-matrix for a fermionic 
contact potential; (b) Relation between the two-particle cor- 
relation function L and the many-particle T-matrix. 

The fermionic correlation functions for the number 
density and current can be quite generally obtained from 
their definitions as thermal averages of the time-order 



operator acting on the number density and current op- 
erators. This is done by expressing the number density 
and current operators in terms of Nambu field operators 
and introducing the two-particle correlation function (14) 
accordingly. One ends up with the following expressions: 

V (a)-~ V r3 r3 fd^dk^±sp 
XnnKQ) - 2^ ^i,e'n»,i'" / (2n)e 132 

e 01 an am J \ ) „i 



xL^ f,„ {k, k'; q) 



(20) 



1 3 fdkdkf 

e,e',e'',e"' ^ ^ 

x(2k'+q)-1^4^f::(fc,fc';g) (21) 

n,n' 

_ 1 f dkdk' 

^^^(^^--(2H^^^^,^/^^^^ 

X (2k + q)(2k' + q)-^ ^ L^'.^k, k'; q) . (22) 

n,n' 

The notation is g = (q, fJ,y), fc = (k, ci;„), fc' = (k',cj„'), 
where t° and are Pauli matrices, and the labels for 
L are specified in Fig. 1(b). The first term on the right- 
hand side of the diagrammatic representation of Fig. 1(b) 
corresponds to the standard BCS contribution to the re- 
sponse functions^^. This term vanishes in the strong- 
coupling limit for temperatures much below the disso- 
ciation temperature of the composite bosons. The sec- 
ond term, on the other hand, corresponds to the broken- 
symmetry BCS-RPA approximation and is required to 
maintain full gauge invariance^^. In the present context 
this term has a nontrivial strong- coupling limit and yields 
the response of the composite bosons condensate. From 
this term one obtains for the correlation functions (20)- 
(22): 

^ De,e{q)Tf/,:(q)De,,,,,,{q) (23) 
',e",e"' 

De,e'iq)Tf/n{q)Ven,^e"{q) (24) 
^ Ve,e'{q)Tllf::{q)Ve'",e'{q) (25) 



Xnn{q) 

Xniiq) 
with 



DiAq)^E-ke"'lj^-^i: 

xQe,e"{k)ge'n^^,{q + k) (26) 

X (2\^ + <i)gi^e'{k)gt" ,t{q + k). (27) 

Recalling that, for a fermionic contact potential, only 
the four elements (16) and (17) of the many-particle T- 
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matrix are nonvanishing, and exploiting the symmetry 
properties of the quantities (26) and (27), one arrives at: 



Xnn{q) = {D{q),D{-q)) 



Xni{q) = {D{q),D{-q)) 



Xjj(9) = (V(g),V(-g)) 



Tiiiq) T,2{q) \ ( D{q) 



T2l{q) T22{q) / V D{-q) 



Tii{q) Ti2{q) \ ( Y{q) 



T2i(g) T22{q) J \ y{-q) 



Tiiiq) Tuiq) 
T2i{q) T22{q) 



V(g) 

Y{-q) 



where D(q) ~ — 4mA(mai?)/(167r) and V(g) ~ 
— qA(mair)/(167r). Using the approximate expressions 

(16) and (17), one obtains finally for the fermionic cor- 
relation functions in the strong-coupling limit: 



Xnn{q) - D^q = 0) [Tn{q) + Tn{-q) + 2T,2{q)] 



= 4 



riB 



TUB {ift,)^ - EB{q)^ 
XnM ^ D{q = 0)V{q) [Tn{q) - T,,{-q)] 
Ub iri^q 



(28) 



ruB (*0„)2 - EBicir 
Xjj(5) - V(rz)V(g) [Tn{q) + Tn{-q) - 2T,2{q)] 



(29) 



_ . Mi3 + q-/''(4;)/_B) 



(30) 



The analytic continuation iO^ —^ui + iS (5 = 0+) to the 
real frequency uj is now required to obtain the retarded 
correlation functions. In this way the bosonic expres- 
sions (6), (12), and (13) are recovered from their fermionic 
counterparts in the strong-coupling limit. There is an ad- 
ditional factor of four, reflecting the fact that the bosonic 
number and current densities are half the fermionic val- 
ues. 

The present mapping has been explicitly established 
for fermions (and composite bosons) in the absence of a 
static external potential. On physical grounds we ex- 
pect the connection wc have established between the 
time-dependent GP equation and the fermionic broken- 
symmetry BCS-RPA approximation to remain valid also 
in the presence of a static confining potential. 

The above forms of the correlation functions obtained 
either from the time-dependent GP equation or from the 
strong-coupling limit of the fermionic broken-symmetry 
BCS-RPA approximation, contain only the contribution 
of the condensate and thus hold at or near zero tem- 
perature. Since the current response of the condensate 
can be only longitudinal, the transverse current corre- 
lation function cannot be obtained in this way since it 
requires inclusion of excitations out of the condensate. 
These excitations are important for determining the be- 
havior of the system at finite temperature and close to 
the superfluid transition temperature. Corrections to the 



time-dependent GP equation, or equivalently additional 
fermionic contributions to the current correlation func- 
tion in the broken-symmetry phase over and above the 
BCS-RPA approximation, arc thus required to recover, 
for instance, the Bogoliubov approximation to the cur- 
rent correlation function for bosons^^. This contains a 
transverse in addition to a longitudinal part. 

In conclusion, we have shown that in the strong- 
coupling limit of the mutual fermionic attraction the re- 
sponse functions obtained at zero temperature from the 
linearized form of the time-dependent GP equation for a 
system of composite bosons are equivalent to those ob- 
tained within the broken-symmetry BCS-RPA approxi- 
mation for the original fermionic system. The present 
mapping extends to the time-dependent case a previous 
result^ that showed the equivalence between the time- 
independent GP equation for composite bosons and the 
strong-coupling limit of the Bogoliubov-de Gennes equa- 
tions for fermions. 

We thank D. Neilson for a critical reading of the 
manuscript. 



^ B. De Marco and D.S. Jin, Science 285, 1703 (1999). 

^ G. Roati, F. Riboli, G. Modugno, and M. Inguscio, Phys. 

Rev. Lett. 89, 150403 (2002). 
^K.M. O'Hara, S.L. Hemmer, M.E. Gehm, S.R. Granade, 

and J.E. Thomas, Science 298, 2179 (2002). 

* For reviews, see E.A. Cornell and C.E. Wieman, Rev. Mod. 
Phys. 74, 875 (2002); D.R. Anglin and W. Ketterle, ibid. 
74, 1131 (2002). 

For reviews, see F. Dalfovo, S. Giorgini, L.P. Pitaevskii, 

and S. Stringari, Rev. Mod. Phys. 71, 463 (1999); A.J. 

Leggett, Rev. Mod. Phys. 73, 307 (2001). 
^ E.P. Gross, II Nuovo Cimento 20, 1766 (1961); J. Math. 

Phys. 4, 195 (1963). 
''L.P. Pitaevskii, Sov. Phys. JETP 13, 451 (1961). 

* C.J. Pethick and H. Smith, Bose-Einstein Condensation in 
Dilute Gases (Cambridge Univ. Press, Cambridge, 2002). 

^P. Fieri and G.C. Strinati, Phys. Rev. Lett. 91, 030401 
(2003). 

P.G. Dc Gcnncs, Superconductivity of Metals and Alloys 
(Benjamin, New York, 1966). 

See, e.g., G.C. Strinati, La Rivista del Nuovo Cimento 
Vol.11, N.12, (1988). 

C.A.R. Sa dc Mclo, M. Randeria, and J. Engelbrecht, Phys. 
Rev. Lett. 71, 3202 (1993); P. Fieri and G.C. Strinati, 
Phys. Rev. B 61, 15370 (2000). 

See, e.g., A.L. Fetter and J.D. Walecka, Quantum Theory 
of Many-Particle Systems (McGraw-Hill, New York, 1971). 
J.R. SchriefFer, Theory of Superconductivity (Benjamin, 
New York, 1964). 

G.C. Strinati, in Electrons and Photons in Solids (Scuola 
Normale Superiore, Pisa, 2001), p. 403. 



4 



